X 



BAR4R-CONF-01/03 
July, 2001 



Measurement of the branching fraction for the decay B° — > D* + D* 

"3 : 

The BABAR Collaboration 
04 ; July 21, 2001 

> 

in 

o 
t> 

o 



Abstract 



Decays of the type B — > D^*'D^*> can be used to provide a measurement of the parameter sin2/3 
of the Unitarity Triangle that is complementary to the measurement derived from the mode B° — > 
J/tpKg. In this document we report a measurement of the branching fraction for the decay B° — > 
D* + D*~ with the BABAR detector. With data corresponding to an integrated luminosity of 20.7 fb~ 
collected at the T(45) resonance during 1999-2000, we have reconstructed 38 candidate signal events 
in the mode B° — > D* + D*~ with an estimated background of 6.2 ± 0.5 events. From these events, 
we determine the branching fraction to be B(B° —> D* + D*~) = (8.0 ± 1.6(stat) ± 1.2(syst)) x 
10 -4 (preliminary) . 
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1 Introduction 

One of the most important goals of the BABAR experiment is to precisely measure the angles 
of the Unitarity Triangle. While the decay B° — > J/ipK® can be used to measure sin2/3, the 
Standard Model predicts that the time-dependent CP violating asymmetries in the decays Q 
B° — > D^*> + D^*> can also be used to measure the same quantity. An independent measurement 
of sin2/3 in these modes would therefore provide a consistency test of CP-violation in the Standard 
Model. 

The vector-vector decay B° — ► D* + D*~ is not, however, a pure CP eigenstate. A sizeable dilu- 
tion of the measured asymmetry may be produced by a non-negligible P-wave CP-odd component. 
The dilution can, in principle, be completely removed by a time-dependent angular analysis of the 
decay products M. 

The rate for the Cabibbo-suppressed decays B — ► D^*'D^*> can be estimated from the measured 
rate of the Cabibbo-favored decays B — > D£'D(*>: 

B(B -► DWSW) « [ ^- ] tan 2 % • B(S -► D^D^), (1) 

where 0c is the Cabibbo angle, and f D (*) and /„(*) are decay constants. From this it follows 
_ u s 

that the B — > D^*>D^*> branching fractions are of the order of 0.1%. Previous measurements of 
branching fractions and upper limits for these modes are summarized in Table |l[ 

Table 1: Summary of branching fraction and upper limit measurements performed by the CLEO 
experiment ||. Upper limits are quoted at the 90% confidence level. 

Decay Branching Fraction (xlO -4 ) 

B u -f D*+D*~ 9.9±%%(stat)±1.2(syst) 
B°^D*+D- < 6.3 

B° -► D+D- < 9.4 



2 The BABAR detector and dataset 

The data used in this analysis were collected with the BABAR detector Q at the PEP-II storage 
ring H located at the Stanford Linear Accelerator Center. The results presented in this paper 
are based on data taken during the 1999-2000 run. This data sample represents an integrated 
luminosity of 23.3 fb~ , with 20.7 fb~ collected on the T(AS) resonance. The total number of BB 
pairs produced in this sample was N B g = 22.7 x 10 6 . 

Charged particles are detected and their momenta measured with the combination of a 40-layer 
drift chamber (DCH) and a five-layer silicon vertex tracker (SVT) embedded in a 1.5 T solenoidal 
magnetic field. Photons are detected by a Csl electromagnetic calorimeter (EMC) that provides 
excellent angular and energy resolutions with a high efficiency for energies above 20MeV. Charged 
particle identification is provided by the specific ionization loss (dE/dx) in the tracking devices and 
by an internally reflecting ring-imaging Cherenkov detector (DIRC) covering the barrel region of 
the detector. 



Table 2: 
K 



o 



IT ' 7T 



D° and D + decay modes and branching fractions R. The branching fraction for 
is included for modes containing a Ki 



-o 



Decay Mode 



Branching Fraction (%) 



TT 



D 



K~TX A 



D° -» K-7T+7T 



3.83 ±0.09 

13.9 ±0.9 

7.49 ±0.31 



L>° -» A"°7r + 7r- 




1.85 ±0.14 


Total D [> Branching 


Fraction 


27.1 


Decay Mode 




Branching Fraction (%) 


£>+ -» ii"°7r + 

L> + -» A'-.K' + 7r + 




9.0 ±0.6 
0.99 ± 0.09 
0.87 ±0.07 


Total -D + Branching 


Fraction 


10.9 



3 Determination of B(B° -> D*+D* ) 

B° mesons are exclusively reconstructed by combining two charged D* candidates reconstructed in 
a number of D* and D decay modes. Events are pre-selected by requiring that there be three or 
more charged tracks and that the normalized second Fox- Wolfram moment || of the event be less 
than 0.6. We also require that the cosine of the angle between the reconstructed B direction and 
the thrust axis of the rest of the event be less than 0.9. 

Charged kaon candidates are required to be inconsistent with the pion hypothesis, as inferred 
from the Cherenkov ring measured by the DIRC and the dE/dx as measured by the SVT and DCH. 
There are two exceptions to this: tighter kaon identification is applied to one of the charged kaons 
in decay D + — > K~K + ir + , and no particle identification requirements are made for the kaon from 
the decay D° -> K~ty + . 

Kg — ► 7r + 7r~ candidates are required to have an invariant mass within 25 MeV/c 2 of the nominal 
K^ mass. The opening angle between the flight direction and the momentum vector of the K® 
candidate is required to be less than 200 mrad, and the transverse flight distance from the primary 
event vertex must be greater than 2 mm. 

Neutral pion candidates are formed from pairs of photons in the EMC with energy above 30 MeV, 
an invariant mass within 20 MeV/c 2 of the nominal 7T° mass, and a summed energy greater than 
200 MeV. A mass-constraint fit is then applied to these 7r° candidates. The 7r° from D* + — > D + tt° 
decays ("soft" n ), however, is required to have an invariant mass within 35 MeV/ c 2 of the nominal 
7T° mass and momentum in the T(4S) frame of 70 < p* < 450MeV/c, with no requirement on the 
summed photon energy. 

The decay modes of the D° and D + used in this analysis were selected by an optimization of 
S 2 /(S + B) based on Monte Carlo simulations, where S is the expected number of signal events and 
B is the expected number of background events. The D° and D + modes used and their branching 
fractions are summarized in Table ||. D° (D + ) meson candidates are required to have an invariant 
mass within 20 MeV/ c 2 of the nominal D° (D + ) mass. 



The D* + mesons are reconstructed in their decays D* 



D°tt + and D* 



D + 7T°. We include 



Table 3: D* decay modes and branching fractions J7|. 



Particle Decay Mode Branching Fraction (%) 

D*+ D* + -> L> u tt + 67.7 ±0.5 

D* + -> D+vr 30.7 ±0.5 

Total Visible Z}* + Branching Fraction 98.4 



for this analysis the decay combinations D* + D*~ decaying to (D° tt + , D° ir~) or (D° ir + , D~ ir°), 
but not (D + 7T°, D~ 7T°) due to the smaller branching fraction and larger expected backgrounds. 
The branching fractions for these modes are summarized in Table |3[ D° and D + candidates are 
subjected to a mass-constraint fit and then combined with soft pion candidates. A vertex fit is 
performed that includes the position of the beam spot to improve the angular resolution of the soft 
pion. 

To select B° candidates with well reconstructed D* and D mesons, we construct a \ 2 that 
includes all measured D* and D masses: 



2 _ (m D - mDp D G V . I m D ~ m D 



XMass — ~ + 



! I I I FT t I I T\ 

PDG 



o«. D 



'rtif 



/Amn» - Amn« \ 2 / AtTT-t^* — A?TT,n* N ' 
I u ^PDG \ I / D U PDG 



\ 



+ 



where the subscript PDG refers to the nominal value, and Am is the D* — D mass difference. For 
<T mD we use values computed for each D candidate, while for a^ m we use fixed values of 0.83 MeV/c 2 
for D* + -> D°7T + and 1.18MeV/c 2 for D* + -* D + ir°. A requirement that xliass < 20 is applied to 
all B° candidates. In events with more than one B° candidate, we chose the candidate with the 

lowest value of XMass • 

A B meson candidate is characterized by two kinematic variables. We use the energy-substituted 
mass, mES> defined as 



m E s = \l Ef eam ~ P* B 2 
and the difference of the B candidate's energy from the beam energy, AE , 

AE = E* B - E* Beam 

where E B {p* B ) are the energy (momentum) of the B candidate in the center-of-mass frame and 
E* Beam is one-half of the center-of-mass energy. The signal region in the AE vs. ities plane is defined 
to be \AE\ < 25MeV and 5.273 < ities < 5.285 GeV/c 2 . The width of this region corresponds to 
approximately ±2.5a in both AE and mss- 

To estimate the contribution from background in the signal region, we define a sideband in the 
Ai? vs. tties plane as 

| AS | < 200 MeV 

5.20 < m ES < 5.26 GeV/c 2 
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and 

50 < \AE\ < 200 MeV 

5.26 <to ES < 5.29GeV/c 2 

We parameterize the shape of the background in the AE vs. toes plane as the product of an 
ARGUS function || in toes an d a first-order polynomial in AE . Based on this parameterization 
we estimate that the ratio of the number of background events in the signal region to the number in 
the sideband region is (1.72 ± 0.10) x 10 -2 . The uncertainty is derived from the observed variation 
of this ratio under alternative assumptions for the background shape in toes and AE . 

Figure p] shows the events in the AE vs. toes plane after all selection criteria have been applied. 
The small box in the figure indicates the signal region defined above, and the sideband is the entire 
plane excluding the region bounded by the larger box outside the signal region. There are a total of 
38 events located in the signal region, with 363 events in the sideband region. The latter, together 
with the effective ratio of areas of the signal region to the sideband region, implies an expected 
number of background events in the signal region of 6.24 ± 0.33(stat) ± 0.36(syst). The quoted 
systematic uncertainty comes from the background shape variation discussed previously. Figure |2| 
shows a projection of the data on to the toes axis after requiring \AE\ < 25 MeV. 
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Figure 1: Distribution of events in the AE vs. toes plane. The small box indicates the signal 
region, while the sideband region is everything outside the larger box. 

We use a detailed Monte Carlo simulation of the BABAR detector to determine the efficiency for 
reconstructing the signal. This, together with the total number of BB pairs produced during data 
collection, allows us to determine a preliminary branching fraction for B° — » D* + D*~ to be 



B(B° -> D* + D*~ 



1.0 ± 1.6(stat) ± 1.2(syst)) x 10" 



The dominant systematic uncertainty in this measurement comes from our level of understand- 
ing of the charged particle tracking efficiency (9.4%). The high charged particle multiplicity in 
this decay mode makes this measurement particularly sensitive to tracking efficiency. Uncertainties 
were assigned on a per track basis for tt, K and slow w, and were added linearly due to large 
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Figure 2: Distribution of events in tties plane with a cut of \AE\ < 25 MeV applied. The curve 
represents a fit to the distribution of the sum of a Gaussian to model the signal and an ARGUS 
function M to model the background shape. 



correlations. The imprecisely known partial-wave content of the B° — ► D* + D*~ final state is an- 
other source of systematic uncertainty (6.6%). This was estimated by calculating the change in 
the reconstruction efficiency for different final angular states in Monte Carlo. Other significant 
systematic uncertainties arise due to the uncertainties on the D* + , D° and D + branching fractions 
(5.6%) and the differences in mass resolutions between Monte Carlo and data (4.1%). The total 
systematic uncertainty from all sources is 14.5%. 

4 Summary 

Using data collected by the BABAR experiment during 1999-2000, we have observed a signal of 
31.8±6.2(stai)±0.4(syst) events in the decay B° — » D* + D*~ . We measure a preliminary branching 
ratio to be 

B(B° -► D* + D*-) = (8.0 ± 1.6(stat) ± 1.2(syst)) x 10" 4 
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